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Abstract. In this paper we study the role of surface plasmon modes in the Casimir effect. First 
we write the Casimir energy as a sum over the modes of a real cavity. We may identify two 
sorts of modes, two evanescent surface plasmon modes and propagative modes. As one of the 
surface plasmon modes becomes propagative for some choice of parameters we adopt an adiabatic 
mode definition where we follow this mode into the propagative sector and count it together with 
the surface plasmon contribution, calling this contribution "plasmonic". The remaining modes 
are propagative cavity modes, which we call "photonic". The Casimir energy contains two main 
contributions, one coming from the plasmonic, the other from the photonic modes. Surprisingly 
we find that the plasmonic contribution to the Casimir energy becomes repulsive for intermediate 
and large mirror separations. Alternatively, we discuss the common surface plasmon defintion, 
which includes only evanescent waves, where this effect is not found. We show that, in contrast 
to an intuitive expectation, for both definitions the Casimir energy is the sum of two very large 
contributions which nearly cancel each other. The contribution of surface plasmons to the Casimir 
energy plays a fundamental role not only at short but also at large distances. 



1. Introduction 

An important prediction of quantum theory is the existence of irreducible fluc- 
tuations of electromagnetic fields even in vacuum, that is in the thermodynamical 
equilibrium state with a zero temperature. These fluctuations have a number of 
observable consequences in microscopic physics for example in atomic physics the 
Van der Waals force between atoms in vacuum. 

Vacuum fluctuations also have observable mechanical effects in macroscopic 
physics and the archetype of these effects is the Casimir force between two mir- 
rors at rest in vacuum. This force was predicted by H. Casimir in 1948 [lj who 
considered two plane parallel perfect reflectors as shwon in Figure [U and found an 
interaction energy i?cas depending only on geometrical parameters, the mirrors 
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Figure 1: The configuration 
studied by H. Casimir 



distance L and surface A> L 2 , and two fundamental constants, the speed of light 
c and Planck constant H 

The signs have been chosen to fit the thermodynamical convention with the minus 
sign of the energy -Ecas corresponding to a binding energy. The Casimir energy for 
perfect mirrors is usually obtained by summing the zero-point energies 4r of the 
cavity eigenmodes, substracting the result for finite and infinite separation, and 
extracting the regular expression ([1]) by inserting a formal high-energy cutoff and 
using the Euler-McLaurin formula [2J. 

The Casimir force was soon observed in different experiments which confirmed 
its existence 0111 [5]. Recent experiments have reached a good precision, in the % 
range, which makes possible an accurate comparison between theoretical predic- 
tions and experimental observations [BJ [7] . 

Casimir considered an ideal configuration with two perfectly reflecting mirrors 
in vacuum. But the experiments are performed with real reflectors, for example 
metallic mirrors which have a perfect reflection only at frequencies below a plasma 
frequency cjp or alternatively for mirror separations much larger than the plasma 
wavelength Ap characetristic for the metal. Accounting for this imperfect reflection 
and its frequency dependence is thus essential for obtaining a reliable theoretical 
expectation of the Casimir force in a real situation. 

The consideration of real mirrors is important not only for the analysis of 
experiments but also from a conceptual point of view. Real mirrors are certainly 
transparent at the limit of high frequencies and this allows one to dispose of the 
divergences associated with the infiniteness of vacuum energy. This point was 
already alluded to in Casimir's papers and an important step in this direction 
was the Lifshitz theory of the Casimir force between two dielectric bulks [El [9]. 
However this idea was fully implemented in theoretical derivations after a quite 
long period. 
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In the limit of small separations L <C Ap, the Casimir effect has another inter- 
pretation establishing a bridge between quantum field theory of vacuum fluctua- 
tions and condensed matter theory of forces between two metallic bulks. It can 
indeed be understood as resulting from the Coulomb interaction between surface 
plasmons, that is the collective electron excitations propagating on the interface 
between each bulk and the intracavity vacuum [1U|, 111] . The corresponding field 
modes are evanescent waves and have an imaginary longitudinal wavevector. At 
short distances, surface plasmon modes are known to dominate the interaction and 
the Casimir energy reduces to [HI [13] 

3 L 

E E v \ = -a— with a = 1.193... (2) 
2 Ap 

Surface plasmons play an important role in many fields of physics. Let us 
only mention as recent examples the surface plasmon assisted enhancement of the 
transmission of light through metallic structures |14[ \TE[ [16] or in the context of 
biomolecular physics, the importance of plasmon fluctuations of a 2-dimensional 
Wigner-like crystal where they can generate an attractive component to the dis- 
persion forces between parallel surfaces |17j . 

Here, we derive and expand in more detail the results previously given in Ref. 
[18] . investigating more closely the influence of surface plasmons on the Casimir 
energy, not only at short but at arbitrary distances. 



2. Casimir energy for real mirrors 

We restrict our attention to the situation of two infinitely large plane mirrors 
at zero temperature so that the only modification of the Casimir formula (pQ) is 
due to the metals finite conductivity. This modification is calculated by evaluating 
the radiation pressure of vacuum fields upon the two mirrors [19] 



E = -EEEf ln(l-^HV^)+c.c. (3) 
e k u 

= A f — V = f°° — 



k 



The energy E is obtained by summing over polarization p=(TE,TM), transverse 
wavevector k = (k x ,k y ) (with z the longitudinal axis of the cavity) and frequency 
cj; k z is the longitudinal wavevector associated with the mode. r£ are the reflection 
amplitudes here supposed to be the same for the two mirrors. 

Imperfectly reflecting mirrors will be described by scattering amplitudes which 
depend on the frequency, wavevector and polarization while obeying general prop- 
erties of stability, high-frequency transparency and causality. The two mirrors 
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form a Fabry- Perot cavity with the consequences well-known in classical or quan- 
tum optics : the energy density of the intracavity field is increased for the res- 
onant frequency components whereas it is decreased for the non resonant ones. 
The Casimir force is but the result of the balance between the radiation pressure 
of the resonant and non resonant modes which push the mirrors respectively to- 
wards the outer and inner sides of the cavity [2UJ. This balance includes not only 
the contributions of ordinary waves propagating freely outside the cavity but also 
that of evanescent waves. These two sectors of ordinary and evanescent waves are 
directly connected by analyticity properties of the scattering amplitudes. 

Expression ([3j) holds for dissipative mirrors as well as for non dissipative ones 
|19j . It tends towards the ideal Casimir formula (pQ) as soon as the mirrors are 
nearly perfect for the modes contributing to the integral. 

The reduction of the Casimir energy (J3j) with respect to the ideal formula (pQ) 
due to the imperfect reflection of mirrors is described by a factor 

m = ^~ (4) 

This factor plays an important role in the discussion of the most precise recent 
experiments. 

3. Mode decomposition of the Casimir energy 

We now recalculate the Casimir energy as a sum over the cavity modes using 
the plasma model for the mirrors dielectric function. 

t \ 1 W P 2vrc 
e(uj) = 1 § , A P = 5 

with co-p the plasma frequency and Ap the plasma wavelength. 

In this case the zeros of the argument of the integrand in ([3]) lie on the real 
axis. In fact, they have to be pushed slightly below this axis by introducing a 
vanishing dissipation parameter in order to avoid any ambiguity in expression ([3]) 
[19] . We may then rewrite ([3]) as a sum over the solutions [u^j m of the equation 
labeled by an integer index m 



rl[iu} 2 e 2ik * L = 1. (6) 

Simple algebraic manipulations exploiting residues theorem and complex integra- 
tion techniques [11] then lead to the Casimir energy expressed as sums over these 
modes 
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The prime in the sum over m signifies as usually that the term m = has to be 
multiplied by 1/2. The sum over the modes is to be understood as a regularized 
quantity as it involves infinite quantities. The upper expression contains as limiting 
cases at large distances the Casimir expression with perfect mirrors and at short 
distances the expression in terms of surface plasmon resonances ([2]). 

The ensemble of modes appearing in Eqn. can be separated into two 
different ensembles. The TM polarization admits propagating cavity modes as 
well as evanescent modes, while the TE polarization allows only propagating cavity 
modes. The first ensemble contains two modes u;+ and u;_ which tend to the usual 
surface plasmon modes at short distances, uj- is always in the evanescent sector, 
while cj_|_ lies either in the evanescent or in the propagating sector depending 
on its parameters as shown in Fig. [2j The second ensemble are propagating 
cavity modes which may have TE or TM polarization. In view of the particular 
behavior of the u;+ mode, we adopt an adiabatic mode definition, where we follow 
it continuously from the evanescent into the propagative sector and attribute the 
whole mode to the surface plasmon modes, which we call "plasmonic modes" 
[18j . This denomination is chosen in order to avoid confusion with the (common) 
definition of surface plasmon modes which defines them as evanescent modes only 
and cuts the u>+ mode into two pieces. Plasmonic modes are thus the cavity modes 
living in the evanescent sector at least for some particular value of their parameters. 
In the same line of reasoning we call all propagative modes minus the propagative 
part of the u;+ mode "photonic" modes. Photonic modes are the cavity modes 
propagating for all cavity length. In the limit L — > oo, all propagative modes 
tend asymptotically to the eigenmodes of the perfect cavity. The frequencies of 
the plasmonic modes w+ and w_ degenerate in the limit L — > oo to the surface 
plasmon frequency uiq for a single interface [21j 



UJ± 



c4 + 2|k| 



c4 + 4|k|' 



The Casimir energy can now be rewritten as 
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(8) 



0) 



plasmonic Contribution (E pl ) photonic Contribution (B ph ) 

Note that both contributions have no physical meaning on their own, i.e. one 
cannot measure them separately. The only observable is the total Casimir energy 
which is the sum of both contributions. We may rewrite the plasmonic contribution 
to the Casimir energy in a more explicit way as follows 



pi 



cA 



'" k ^(a; + [k,L]+^[k,L]-2^[k]) 



(2vr) 2 2 



(10) 
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Figure 2: A representa- 
tion of the dispersion rela- 
tion for the TM-modes with 
the plasma model for L/ Ap = 
1.5. The dashed curves repre- 
sent the perfect cavity modes. 
The solid grey curves rep- 
resent the plasmonic modes 
and the black curves the 
photonic modes. The lo + - 
mode changes its nature when 
crossing the light- dispersion 
curve. 



The basic idea of the explicit calculation of the plasmonic contribution to the 
Casimir energy resides in the fact that the frequencies functions u>i, i = 0, ± 
are solutions of simple equations. We rewrite Eq. (110p in terms of dimensionless 
variables 

180 f°° -v 

E p \ = ^pi-Ebas, Vp\ = 3- / ^Cik tti[k]dk (11) 

71 ^° i 

n = ujL, n P = uj P L, |k|L = fc, z = k 2 -n 2 (12) 

with c+ = c_ = 1, Co = — 2. is the dimensionless surface plasmon frequency for 
a single mirror and we have introduced the corrective factor 7] p \ for the plasmonic 
contribution to the Casimir energy. Note that 

180 f°° 

r] t = T knMdk, i = ±,0 (13) 

Jo 

is divergent despite the convergence of the whole expression given in Eq. (jlOf) . 
Without giving any details let us just mention that in order to perform the explicit 
calculation we need to introduce a regularizing factor. Such a modification is 
mathematically mathematical convenient and does not affect the final result. 

It can be shown that the dimensionless frequencies £li[k] can be formally ob- 
tained as 

k 2 = fi(z) => tti[k] = y/tf-fr 1 ^], i = 0, ± (14) 
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where 



f+(z) = z+ ^= = z + gl[ z ] (15a) 



/- (z)=z + fi^E = z + g 2 _ [z] (15b) 

y/z + Jz + nl coth[^] 



SllJz 

fo(z) = z + p r = z + 5o 2 [z] (15c) 



' : + • <>[> 



Let us also define 



v 2 = z° 



-/r 1 [o] = ^ 2 [o] (i6) 

With the change of variable k 2 = z + g 2 (z) and after some rearrangements Eqn. 
(|11|) can be rewritten as 



180 



/ ^2cigi(z)dz + ^g + (z)dz--y% 

JO i J-z\ 



(17) 



where we have exploited the fact that 

z° + ^ z°_ = z° = 0, /r 1 [oo] = oo 
fi(-z?) = 9i (-z?) = ^Jz] and 9i (oo) = ^| 

The corrective factor 77 p i has a well defined structure: it is indeed decomposed 
into an integral over the positive real z-axis plus an integral over an interval of 
the negative z-axis plus a constant depending only on S7p. Moreover only g + , 
is involved in the last integral and in the constant. This particular structure 
can be traced back to the properties of the plasmonic modes. The positive z- 
value domain coincides with the evanescent sector while the negative one describes 
the propagative sector. While the plasmonic mode cu- and the surface plasmon 
frequency for a single mirror ujq are totally contained in the evanescent sector, 
the plasmonic mode u;+ lives in both sectors. Therefore go and g- describing the 
properties of oj- and ujq are contained only in the first integral while g+ has to be 
evaluated in a wider range of z- values which includes an interval in the propagative 
sector. The second integral in Eq. (|17p is thus basically the propagative part 
contribution of the plasmonic mode w+. 

Fig. [3] shows the numerical evaluation of Eqn. (|17p for r] p \ as function of L/Ap 
for two different distance intervals. The left graphic illustrates the short distance 
behavior In the limit f2p <C 1 the corrective factor r] p \ can be approximated by the 
first integral of Eq. (fT7l) . leading to 
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Figure 3: The normalized plasmonic mode contribution to the Casimir force as function 
of L/ Ap for two different distance intervals. 



180 ftp 



rexj 

J \y/l + e"v^ + VI - e -v^ - 2 



3 L 
2°^ 



(18) 



2vT 3 V2. 

We recover here the result of Eqn. ([2]). The right graphic in Fig. [3] shows the 
plasmonic mode contribution rj p \ at large distances. Surprisingly, it changes its sign 
for ~ 0.08 and its slope and diverges for L > Ap [18]. In the large distance 
limit ftp > 1 we find the following asymptotic behavior 



-rvftp, 



29.752 



(19) 



The contribution of the plasmonic modes to the Casimir energy becomes thus 
repulsive for intermediate and large distances. This result is based on an adiabatic 
definition of the surface plasmon modes, where we follow the w+ mode even when 
it crosses the barrier u = ck and becomes propagative. 

Let us now compare this result to the common definition of surface plasmon 
modes which includes only evanescent waves and cuts the w+ mode at lo = ck as 
for example done by Bordag recently [22]. This leads to 



where 



180 

TV 3 



kp 



k(n+[k] - n [k])dk 



180 

7T 3 



fc(n_[fc] - n [k\)dk (20) 



k P =g+(0) = ft P /Vl + ftp/2 (21) 

is associated with the value of |k| for which this modes crosses the light cone. rj ev 
can be evaluated numerically [22], or using the same method developed here 



Vc 



180 ( f°° ■>, 2 

■^3 ( / ^ Ci 9i( z ) dz -J p 9o(z)dz- -(k P -n 3 [k P }) 



where we have exploited the following relations 



- Zn 



fo^kp] = 4- n 2 [k P ] go(-z^) = n [k P ] 



(22) 



(23) 
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Eq. (|22p shows the short distance asymptotic behavior of the total Casimir energy 
given in (J2j) as it was reported in [121 E3J- It does not change its sign and, in 
agreement with [22], at long distances goes as 

r/ ev = /3 ev \/^P wi t h Aw = 1.62399... (24) 

With this definition we therefore naturally find the result that the contribution of 
the evanescent modes to the Casimir energy is always attractive and reproduces 
well the short distance behavior of the Casimir energy. 

4. Conclusion 

As a concluding remark we would like to stress that, no matter how we attribute 
the propagative part of the u;+ mode, whether to the surface plasmon modes in an 
adiabatic definition (plasmonic modes) or to the propagating modes, the influence 
of surface plasmons is very important at all distances. The Casimir energy is 
the detailed balance of two very large contributions of opposite sign which nearly 
cancel each other, but not quite, the difference being a small Casimir energy, much 
smaller than each of both contributions. 

It might be interesting to investigate if a change in the photon-plasmon cou- 
pling could somehow influence this detailed balance and therefore the value or even 
the sign of the Casimir force. A different coupling could be obtained by using for 
example nanostructured surfaces. For such an analysis the adiabatic mode defini- 
tion should be well suited, because for a different coupling the mode will change as 
a whole and in following continuously the mode one could trace back the changes 
introduced through the structured surface. 
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